Spatial search on a honeycomb network 



The spatial search problem consists in minimizing the number of steps required to find a 
given site in a network, under the restriction that only oracle queries or translations to 
neighboring sites are allowed. We propose a quantum algorithm for the spatial search 
problem on a honeycomb lattice with A'^ sites and torus-like boundary conditions. The 
search algorithm is based on a modified quantum walk on an hexagonal lattice and the 
general framework proposed by Ambainis, Kempe and Rivosh HAmbainis et al. 2005 |l is 
employed to show that the time complexity of this quantum search algorithm is 
O(VMogjV). 
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1. Introduction 

Quantum Walks (QW) are useful tools to generate new quantum algorithms (jAmbainis 20041 
IShenvi et al. 200^ lAmbainis et al. 2005^ . For example, the optimal algorithm for solv- 
ing the element distinctness problem, which aims to determine whether a set has re- 
peated elements or not, is based on QWs (jAmbainis 2004p . An optimal search algorithm 
equivalent to the celebrated Grover's algorithm (jGrover 1996p . uses a modified QW 
on an n-dimensional hypercube to find an element among N sites after 0{'\fN) steps 
([Shenvi et al. 2003)) . Although the QW is a unitary (i.e. invertible) process, it is often in- 
troduced as the quantum analog of a random walk or, more generally, of a Markov process. 
There are two versions of QWs: discrete-time (jAharonov et al. 1993[) and continuous- 
time (jFahri and Gutmann 1998)) walks. The first one uses an auxiliary Hilbert space, 
which plays the role of a quantum "coin" whose states determine the directions of mo- 
tion. Even though both types of QW's have similar dynamics, they are not equivalent. 
For instance, the optimal algorithm for spatial search in two-dimensional grids using the 
continuous-time version has no advantage over the classical algorithm in terms of time 
complexity (jChilds and Goldstone 2004]) . while the algorithm based on the discrete-time 
version has an almost quadratic improvement (jTulsi 2008^ . 

Grover's algorithm applies to non-ordered databases, where there is no notion of dis- 
tance between two elements. However, when storing information in physical memory, 
a given item is stored at a specific location. This poses an interesting alternative ver- 
sion of searching, called spatial search, as the problem of finding a marked location 
in a rigid structure using only local operations: in one time step one can either query 
an oracle for the given site or move to a neighbouring site. BeniofF (jBenioff 2002^ ad- 
dressed this problem on a two-dimensional square lattice with N points. He was the first 
to point out that a straightforward application of Grover's algorithm with the spatial 
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search constrain requires fl(N) steps with no improvement over classical algorithms in 
terms of time complexity. Aaronson and Ambainis (jAaronson et al. 2003P have developed 
a quantum algorithm for this problem with time complexity 0{^/N \og^ N). Ambainis, 
Kempe and Rivosh (AKR) (jAmbainis et al. 2005^ have proposed a QW-based algorithm 
which improves the time complexity to 0{VN log N). Recently, Tulsi (|Tulsi 2008)) has 
proposed an improved version of the AKR spatial search algorithm for two-dimensional 
square lattices, with a time complexity of 0{^/N logN). It is an open problem whether 
the lower bound n(-\/iV) can be achieved for the spatial search on two-dimensional lat- 
tices (jBennet et al. 1997| . AKR have proposed a generalized framework for QW-based 
algorithms on lattices of arbitrary structure, in which the time-complexity of the algo- 
rithm may be obtained from the eigenvalue spectrum of the QW evolution operator. 
Following AKR, we shall refer to this formalism as the abstract search framework. 

In this paper, we provide a new QW-based algorithm which solves the spatial search 
problem in a hexagonal (honeycomb) network in 0{\/N log N) steps. The time com- 
plexity is analyzed using the abstract search framework just discussed. The hexagonal 
network has received attention from condensed matter physicists for many years, due 
to its role in the band theory of graphite () Wallace 1947[) . More recently, the develop- 
ment of graphenes (two-dimensional hexagonal arrays of Carbon atoms) and its possible 
uses in quantum computation (|Van den Nest et al. 2006P have renewed the interest on 
these networks (jGeim et al. 2007|) . The paper is organized as follows. In Section [5] we 
discuss the implementation of a quantum walk on a periodic hexagonal network and ob- 
tain the evolution operator in the Fourier-transformed space. In Section [3] we summarize 
the abstract search framework and use it to evaluate the time complexity of the search 
algorithm on a hexagonal lattice. In Section 2] we present our conclusions. 

2. QW on the hexagonal network 

The Hilbert space of a QW, H — He <E)'Hp is composed of a coin. He, and a position 
subspace. Hp. The evolution operator is of the form U = S ■ {C(E)I) where C is a unitary 
operation in He, I is the identity in Hp and S, a shift operation in H, performs a condi- 
tional one-step displacement as determined by the current coin state. The main challenge 
to obtain the time complexity of a QW-based algorithm on a honeycomb lattice is the 
calculation of the spectral decomposition of the evolution operator U of the underlying 
QW. The abstract search framework is based on a modified evolution operator U' = S-C, 
obtained from the standard quantum walk operator U by replacing the coin operation 
C with a new unitary operation C" which is not restricted to He and acts differently 
on the searched vertex. Ambainis and coworkers have shown that the time complexity 
of the spatial search algorithm can be obtained from the spectral decomposition of the 
evolution operator U of the unmodified QW (jAmbainis et al. 2005| . which is usually 
simpler than that of U' . 

In regular networks, the use of the Fourier transform on the spatial coordinates con- 
siderably simplifies the expressions for the eigenvalues and eigenvectors. It is known that 
a Bravais lattice has an associated reciprocal lattice (|Kittel 1995P and this provides a 
systematic way for obtaining the Fourier transform. The honeycomb network is not a 
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Figure 1. Elementary vectors for the honeycomb network. The white sites form a 
lattice and black sites form the associated basis. 

Bravais lattice, but this can be circumvented by splitting the vertices into two sets with 
N/2 sites each (the lattice and basis sets) and encoding the which-set information on an 
auxiliary one-qubit state. In Fig. [H we distinguish between the N/2 lattice sites (gray) 
and the N/2 basis sites (black) using a color code. 

Let us consider the distance between two adjacent sites of the hexagonal network as 
the unit distance. Then, the vectors ai and a.2 which connect two neighboring lattice 
sites (see Fig.[T]) have norm VS and span an angle of 60°. The unit vector b which locates 
the basis site adjacent to a given lattice site is given by b = ■|(ai + a2). An arbitrary 
lattice point may be addressed by a vector with integer components 



and each lattice point has an associated basis point at r + b. Assume periodicity in both 
directions (von Karmann boundary conditions), so that ni,n2 G [0, m — 1] . For simplicity, 
we consider a number of sites N such that N ~ 2m^, for some integer m. Thus, for an 
A^-element network, we have N/2 kets |rti, 71,2) spanning the position subspace associated 
to the lattice. The N/2 basis sites are accounted for by introducing an auxiliary qubit, 
{|0), |1)}, which is zero for a lattice site and 1 for a basis site. Thus |0;ni,ri2) = |0) (8) 
|ni,ri2) indicates a state associated to a lattice site and |l;rti,7i2), the state associated 
to the corresponding basis site. The TV-dimensional lattice subspace. Hp, is spanned by 
kets {|s; ni, ni)} with s = 0, 1. 

At a given site there are three possible directions of motion. We label each of them 
with an integer index j = 0,1,2 so that the direction of motion is encoded in a three- 
dimensional "coin" subspace. He, spanned by {|0), |1), |2)}. The full 3A^-dimensional 
Hilbert space is H — He ®Hp and the basis states |j;s; 711,712) form an orthonormal 
set. In this basis, a generic state l^*) e "H is expressed as 



where the complex coefficients aj,„j.„2 {bj,ni.n2) ^-re the lattice (basis) components and 
the normalization condition (^'|^') = 1 is assumed. A step in any direction from a lattice 
(basis) point leads to a basis (lattice) point, according to the propagation rule 



r = Tliai -I- 77232 



(1) 





|j;s;77i,772) \j; s ® 1; m ~ {-lyaj, 772 - {-iy(3j) 



(3) 
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where © is the binary sum and Vj — {aj,l3j) are the directional vectors 

vo = (0,0), VI = (1,0), and V2 = (0,l). (4) 
This conditional displacement is implemented with a shift operator, 

S=T. s © 1. ft - (-1)' vj)0-, s, fi| (5) 

j, s, n 

where we have introduced the shorthand notation n for (ni,n2) and the sum modulo m 
is understood for these components. The evolution operator of a quantum walk on the 
hexagonal network is then 

U = S-{G3<» Ip) (6) 

where Ip is the identity in "Hp. The three-dimensional Grover operation G3 acts in T-Lc 
and, in the representation stated above, is given by 

2 \ 

(7) 

After t iterations, an initial state |5'o) evolves to \^t) = J7*|\l/o)- Note that ?7 is a real 
operator, as required by the abstract search formalism (jAmbainis et al. 2005p . 

For single-step displacements, the spatial part of the evolution operator is diagonal 
in the Fourier representation, so let us now consider the Fourier transform in T-Lp. The 
reciprocal lattice (jKittel 1995^ to the one defined by the vectors {ai,a2} is formed by 
vectors {gi, §2}, which satisfy 

gi-ai = g2-a2 = 27r/m, 

gi • a2 = g2 • ai = 0. (8) 

A point of the reciprocal lattice is located through a vector k = fcigi -I- fc2g2 for integers 
ki,k2 G [0, m — 1]. We shall use the short-hand notation k for the two-component vector 

(fel,fc2). 

The coin components of |^) play no essential role in what follows, so let us for the 
moment omit the coin dependence. Then, a state j^*) can be expressed either in the 
position representation or in the wavenumber representation as 

|*> = ^ (a^ilO; n) + bf^l; n)) = ^ (/j|0; k) + g^\l; k)) . (9) 

n k 

The N states |s; k) are related to the position representation by the Fourier transform 

i^;'^) - \/f E^^^'^'H^^ft)- (10) 



s;n = 



lE^^'-^'k). (11) 



N 

k 



These states satisfy (s,k|s',n) = y e*'' '' (5s,s', so Fourier transformed kets of lattice 
(basis) states are orthogonal to basis (lattice) kets. 
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Taking into account the coin dependence and using the above relations, the action of 
the shift operator, eq. ([5]), on k-space is 

5|j;s;k) -c.-(-i)'''^-Mj;^*®l;k), (12) 

where w = exp(27rj/m) and the directional vectors Vj have been defined in eq. (jlj. Notice 
that 5* is diagonal in /c-space and connects lattice points with basis points as expected. 
This fact effectively reduces the problem to a six-dimensional subspace >Ck spanned by the 
kets {|j; s)}. Since k takes N/2 values, the Hilbert space is now decomposed in this sub- 
space and the one spanned by the |k) states, with a dimensional count 6 x N/2 = 3A^. In 
this six-dimensional subspace, in the representation stated above, the reduced evolution 
operator J/k has the explicit form 
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(13) 



Its characteristic polynomial factors as 

P(A) = (A - 1)(A + 1)(A4 - 2 cos(20fc) A^ + 1), (14) 
where the angle 9k G [0, ^] is defined by 

cos(20fc) = - I cosfci + cosfc2 -I- cos (fci — ^2) I — -r, (15) 
9 V / 3 

and ki = 27rfc,j/m for i = 1, 2. The six eigenvalues of f7k are ±1 and ±6^*^*=. 



3. Time complexity of the search algorithm 

The abstract search formalism described in (|Ambainis et al. 2005^ provides a way to im- 
plement a spatial search algorithm on a network where a QW has been properly defined. 
A convenient summary of the abstract search formalism can be found in Ref. ()Tulsi 2008p . 

Assume that the search is for a single site, r = ro, in a periodic hexagonal (honey- 
comb) network with TV sites. The effective target state in Ti. is \t) = \u) ® \tq) where 
|w) = "Y^j g |j, s) is the uniform superposition in £k- 

The generalized search algorithm iterates the unitary operator 

U'^U-Rt, (16) 

where U is the unperturbed quantum walk operator defined in eq. ([6|) and Rt = I3N — 
2\t){t\. In the introduction, we mentioned that a generalized search is implemented with 
a modified quantum walk operator of the form U' — S ■ C , where C is a unitary coin 
operation that acts differently on the searched site, i.e. C' = C ® [Ip — |ro)(ro|) + Ci® 
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|ro)(ro|. Both forms for U' are equivalent, provided the Grover coin C = G3 is used and 
the usual choice of Ci — ~Ic is made for the coin operation on a searched site. 



The initial state for the algorithm is the uniform superposition in "H, 



\'^o)^\u)®\up)^-=Y,\r,s;h), (17) 



where \up) = y |n) is the uniform superposition in position space. Except for a 

phase shift, the operator Rt implements a reflection about the effective target \t) and 
a single application of Rt on the uniform superposition "marks" the searched state by 
changing its relative phase, in a similar form as in Grover's search algorithm (jGrover 1996[) . 

As mentioned previously, Ambainis et al. prove the remarkable result that the time 
complexity of the abstract search algorithm depends on the eigenproblem of U alone (jAmbainis et al. 2005|) 
They show that, after T = 0{l/a) iterations of Ua, the initial state evolves to a final 
state 1$/) = J/Jl^o) which has an increased overlap with the effective searched 

state \t). Detailed expressions for the dependence of a and {^f\t) on the eigenvalues 
and eigenvectors of U are given below. The unperturbed operator U must satisfy two 
conditions: (i) U must be a real operator and (ii) the uniform superposition state |$o) 
must be a non-degenerate eigenstate of U with eigenvalue 1. Both conditions are met by 
the quantum walk operator U defined in eq. since G3 is real and (G3 (E) -^2)!") = \u). 

We follow the notation of Ref. ()Tulsi 2008^ to describe the eigenproblem for U. The 
eigenvectors associated with the —1 eigenvalue, which may be M-degenerate, are labeled 
as for i = 1 . . . M. Let |$^) indicate the eigenvectors associated to all other eigen- 
values distinct from ±1. The eigenvectors may be chosen so that the amplitudes on \t) 
on the proper basis of U are real. Then, the effective target state may be expanded with 
real coefficients as 



where the index £ runs over all pairs of conjugate eigenvectors with eigenvalues distinct 
from ±1. These amplitudes {ao,ae,ai), together with the angles 9k defined by eq. p5|) . 
determine the time complexity of the abstract search algorithm (jAmbainis et al. 2005| 
ITulsi 2008| . The rotation angle towards the searched element, which results from a single 
application of U\ is 



M 




(18) 




(19) 



After T = 7r/2a iterations, the overlap with the searched state is 




(20) 



In both expressions, the sums run over the eigenvalues distinct from ±1. 
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The (unnormalized) eigenvectors \v^) associated with the eigenvalues ±1 are 



1 



k-2 



icj'^^ (1 - 
w''! - 1 



(21) 



except for k\ — k2 — 0. Note that the projection of the effective target state |t) on £k 
is the uniform state \u) and {u\v^) = 0, unless k\ = ki = 0. In this degenerate case, 
the eigenvalues are ±1 and |u) itself is an eigenvector of C/k with eigenvalue +1. All the 
other eigenvectors are orthogonal to \u) so, for all k, 



k.O 



and 



= 



(22) 



so that flo = y^/iV and the terms corresponding to the eigenvalue —1 do not contribute 
in eq. (fT8|) . Let us indicate the eigenvectors associated to the other eigenvalues ±6"*='^* 



as I ± j^, 



Then, eq. ([T8| for the effective target state reduces to 



K> - + \/|e K (i + -t') + 1 + -^-'^''>) + 

k#0 



with the real amplitudes 



«k(|-4'''^) + |-^^''"^))]®|k), (23) 



1 + cos ki + cos fc2 



'1± 



3 cos 9k 



(24) 



Even though analytical expressions for all the eigenvectors of U are unknown, knowledge 



of the coefficients allows us to evaluate the time complexity of the search algorithm. 



For the quantum walk on a honeycomb, eq. (|19p leads to 
1 



= O 



St 

k#0 



— cos 9k 1 + cos 6*^; 



(25) 



Let us concentrate on the A^-dependence, for N ^ 1, of the argument A{N) of the above 
square root. Using eq. ()24p . after some manipulation, we obtain 



A 



-E- 

6kVo 



cos fci + COS A:2 
sin^ 6k 



N 



48 (tt - ey 



dk2dki 



- 4 + cos fci + cos ^2 



sin^ 0fc 



(26) 



where we have used sin^ 6k — i 



^cos fci + cos ^2 + cos(fci — ^2)^ and approximated 



the sum by an integral in the usual form, X^k^^c 



^ f (7F^ / /r ' dkdh with £ = 27rVW- 
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For 3> 1 (or e <C 1), the A^-dependencc of A is 

^ ' 32 J, J, kl + kl - k,k2 He' 



So, l/a = 0(V^/VloglV) iterations of U' are required to reach the final state 

Using eq. ([20| . we obtain that the inverse of the overlap between the final state and 
the target \t) is 



'f/)! = O Ut V (4) + K) cot^ (^^/4) ^ O {B{N)) . (27) 




Using eq. dM]) and for iV > 1, the iV-dependence of B{N) is 

E cot' (^'^■/4) =^ ^(^^ ' *2*i cot^ (^±^^\ogN, (28) 

where the divergence comes, as before, from the sin~^ 6k term. Then 

1 



= 0(logiV). (29) 

i(ti$/>r 

The analysis of the time complexity of the algorithm is as follows. After 1/a = 
0(y/N\ogN) iterations of U' , the algorithm reaches the final state |$/) with probabil- 
ity p = . The method known as amplitude amplification ("Br assard et al. 2002^ 
states that if there is an unitary operator U' such that the probability of measuring a 
marked state upon measuring C/'*|$o) is p > 0, then there is a quantum procedure that 
finds the marked state with certainty using C>(l/y^) applications of C/'*. That procedure 
uses the inversion about the mean, which can be implemented in 0{^fN) steps. This 
leads to an overall complexity of O^sfN logiV) to find the marked state in the honey- 
comb lattice. This is the same complexity of the AKR spatial-search algorithm on the 
cartesian grid of a torus (jAmbainis et al. 2005p . where each site has four neighboring 
sites and the N sites form a lattice. 

In a remarkable paper, A. Tulsi (jTulsi 2008^ has described a method to improve even 
further the probability of finding the marked vertex. Let us consider a quantum circuit 
that implements operator Rt followed by U as defined in eq. (|16p . Tulsi introduced an 
extra qubit and defined a new one-step evolution operator as described in the circuit of 
Fig. [21 where — Z is the negative of Pauli Z operator and 

X.-{ (30) 
\ — smd cosd / 



where 5 must assume the value 1 / yloglV. It is straightforward to show that Tulsi's proce- 
dure increases the overlap between the final state and the target, such that |(i|<I>j)| — 0(1). 
Consequently, the overall time complexity of the search algorithm in the honeycomb lat- 
tice may be improved to 0{\/ N\ogN)^ as in the AKR case, with Tulsi's modification. It is 
not necessary to use the amplitude amplification method in this case. We have performed 
an independent numerical simulation which agrees with this analytical calculation. 
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Figure 2. Tulsi's circuit diagram for the one-step evolution operator of the 
quantum walk search algorithm. 

4. Conclusions 

Hexagonal networks (honeycombs) are the underlying representation of a carbon struc- 
ture called graphene, which has been attracting special attention over the last years, 
especially for its potential applications in nanotechnology. In this paper, a new quan- 
tum algorithm for spatial search in a honeycomb with periodic boundary conditions 
is discussed. The protocol is based on a quantum walk in the honeycomb. We obtain 
the expression for the evolution operator in the Fourier representation and solve its 
eigenvalue problem. Then, the abstract search formalism developed by Ambainis et al. 
(jAmbainis et al. 2005P is used to obtain the complexity of the algorithm from the par- 
tially known spectral decomposition of the evolution operator. Our results have been 
verified by numerical simulations. 

The search algorithm on the honeycomb has an overall time complexity of 0{^/N log N) 
by using the amplitude amplification procedure. A better improvement, to 0{^/ N log N), 
can be obtained by using Tulsi's technique. Surprisingly, this is the same complexity found 
for the quantum search on the square grid after Tulsi's improvement. Both the hexagonal 
grid and the square grid are regular graphs which cover the plane, although the former 
has degree 3 and the latter has degree 4. The fact that the complexity of the search 
algorithm is the same in both cases suggests that the number of connections of each 
node is not affecting the complexity of the abstract spatial search algorithm. 

Several open questions remain. One of them is whether the abstract search algorithm 
has the same complexity when applied to graphs of general degrees. The triangular 
network, for instance, has degree 6 and also covers the plane. It would be interesting to 
investigate the behavior of the algorithm on this topology. One may also inquire about 
how robust the search algorithm is when there are some missing nodes. Finally, we point 
out that an optimal spatial search algorithm 0{\/N) for the case of a two-dimensional 
network covering the plane has not yet been found. 
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